Gravitational radiation from collapsing magnetized dust 



Hajime Sotani p Shijun Yoshida and Kostas D. KokkotasH 

^Department of Physics, Aristotle University of Thessaloniki, Thessaloniki 54124, Greece 
^Science and Engineering, Waseda University, Okubo, Shmjuku, Tokyo 169-8555, Japan 

(Dated: February 7, 2008) 

In this article we study the influence of magnetic fields on the axial gravitational waves emitted 
during the collapse of a homogeneous dust sphere. We found that while the energy emitted depends 
weakly on the initial matter perturbations it has strong dependence on the strength and the distri- 
bution of the magnetic field perturbations. The gravitational wave output of such a collapse can be 
up to an order of magnitude larger or smaller calling for detailed numerical 3D studies of collapsing 
magnetized configurations. 
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^ ; I. INTRODUCTION 

(D ...... . 

^ , The direct detection of gravitational waves will be a significant breakthrough both for fundamental physics and 

astrophysics. With the information collected from gravitational waveforms, one will be able to make validation of 
general relativity, collection of astronomical data, and examine the nature of matter in supranuclear densities. Among 
the most important applications of the gravitational wave observations is the asteroseismology in which seismic infor- 
I ' mation for determining stellar structures is obtained by gravitational waveforms [11, [3, [^, 0, 0,0, 0| ■ Currently, several 
J> ground-based acoustic and laser interferometric detectors for gravitational waves like LIGO, TAMA300, GEO600, 
T— I ■ and VIRGO are in operation but there is not yet any direct detection of gravitational waves [1] . In addition to the 
0^ ground-based detectors, there is a project to launch a Laser Interferometric Space Antenna (LISA), which is planned 
■ to be launched by the end of next decade and to operate for five years @ . 

, The nonspherical stellar collapse is one of the potential sources of gravitational waves both for the ground-based 
and for space detectors. The ground-based interferometers target the formation of stellar mass black holes or neutron 
\ stars because they are most sensitive to gravitational wave frequencies in the range 10 - 1000 Hz. Space detectors 
with sensitivities ranging from 10"'* up to 10~^ Hz might detect signals from the creation of supermassive black-holes 
[Tol . [Til . [13I . [3, [13 • There are two approaches to calculate the energy radiated away as gravitational radiation 
r ' during stellar collapse; the first is via direct numerical integration of the exact Einstein and matter equations, which 
^jPj' form a coupled nonlinear system of partial differential equations, and the second is by making use of the linear 
• • , perturbation analysis. During the last decade numerical relativity made remarkable advances and many complicated 
matter and spacetime configurations can be treated with high degree of confidence 

In spite of the recent great progress in numerical relativity, the accurate extraction of gravitational waveforms has 
' been proved a quite difficult task. The reason is that the weak gravitational waves emitted during the stellar collapse 
' sometimes have amplitudes similar to those of unphysical noises due to gauge modes and/or to numerical errors. Linear 
perturbation theory is extremely efficient for such processes since the nonradial perturbations which are responsible for 
the emission of the weak gravitational waves are separated from the nonradiative symmetric background. For example, 
for spherically symmetric backgrounds the physical quantities described by the perturbations can be expanded in 
terms of tensor spherical harmonics and this separation of variables simplifies considerably the study of the problem. 
Actually, the master equations for the perturbations are reduced to simple set of coupled linear partial differential 
equations, which can be evolved with extremely high accuracy. On the other hand via this procedure one might miss 
certain nonlinear phenomena that take place at the very last stages of collapse. 

Linear perturbation analysis has been used in calculations of the gravitational wave emission from the stellar 
collapse to a black hole [12, [H, [13, [1H|- Cunningham, Price, and Moncrief derived the perturbation equations on 
the Oppenheimer-Snyder solution, which describes collapse of homogeneous dust [1^ and calculated gravitational 
radiation emitted during the collapse to a black hole [22| . By using the gauge invariant perturbation formalism on the 
spherically symmetric spacetime formulated by Gerlach and Sengupta [27[ , Seidel and co-workers [23,] investigated the 
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gravitational waves from the stellar collapse in which a neutron star is born. The gravitational waves from collapse of 
an inhomogeneous dust, which can be described b y n onradial perturbations of the Lemaitre-Tolman-Bondi solution 
[28l |. are computed by Iguch, Nakao, and Harada [23] ■ Harada, Iguchi, and Shibata [1^1 calculated the axial parity 
gravitational waves emitted from collapse of a supermassive star to a black hole by employing the covariant gauge- 
invariant formalism for nonradial perturbations on spherically symmetric spacetime and the coordinate-independent 
matching conditions at stellar surface, devised by Gundlach and Martin-Garcfa (29j. 

Although, as mentioned above, there have been many investigations of gravitational radiation from stellar collapse 
with the linear perturbation analysis, effects of magnetic fields on the gravitational radiation have not been taken into 
account. While Cunningham et al. dealt with electromagnetic perturbations of the Oppenheimer-Snyder solution 
[2^ . they omitted the conduction of the fluid. In other words, Cunningham et al. did not consider the direct coupling 
between fluid and magnetic field. However, it has been recently realized the importance of effects of magnetic fields 
on the evolution of compact objects again due to the advent of high-performance instruments like satellite-borne 
detectors. One of the most impressive examples is the discovery of magnetars, which are neutron stars whose strength 
of magnetic fields is estimated about 10^^ G. The magnetar model, in which observed activities are powered by decay 
of the strong magnetic fields of magnetars, successfully explains activity of soft gamma-ray repeaters (SGRs). All 
the four of the known SGRs have rotation periods of '--^ 5 — 8 s, and the three of them have large period derivatives 
of ^ lO^^'' s s^^, which infer the existence of magnetic fields of B ^ (5^8) xlO^'' G [lOl- Since there is an ultra 
strong magnetic field in some neutron stars, which will be born by stellar collapse, it is natural to take into account 
its effect on stellar collapse. Even if the initial magnetic field is weak, it is conceivable that the magnetic fields of 
the collapsing object are, due to the magnetic flux conservation, amplifled during the collapse and will most probably 
affect the gravitational waves emitted. 

Another example for showing the importance of magnetic fields in the evolution of compact objects is related to 
gamma-ray bursts (GRBs). The short-duration GRBs could result from hypergiant flares of magnetars associated 
with the SGRs [31| or magnetized hypermassive neutron star collapse |32j. For long-duration GRBs, strong magnetic 
flelds provide the excitation energy on the required time scale and drive collimated GRB outflows in the form of 
relativistic jets [SS*]. Additionally, the so-called hypermassive neutron star, which can be formed after the merger of 
binary neutron stars and can be in an equilibrium state due to differential rotation, could lead to delayed collapse 
by the magnetic braking and viscosity, even if the initial magnetic field and the viscosity are very weak [s^ . Notice 
that even if they are weak initially, the magnetic fields can be amplified to the required strength by the winding-up 
of weak magnetic field due to differential rotation |35|, |36| . 

All the examples mentioned above indeed suggest that magnetic fields play an important role in the stellar collapse. 
In this paper, we therefore consider gravitational radiation from collapse of weakly magnetized dust spheres to 
explore the effects of magnetic flelds on the gravitational waves from the stellar collapse to a black hole. For the 
sake of simplicity, we consider, as the background spacetime, the Oppenheimer-Snyder solution, which describes 
homogeneous dust collapse. The weak magnetic flelds of the dust spheres are treated as small perturbations around 
the Oppenheimer-Snyder solution. We therefore regard both the gravitational waves and the magnetic flelds as 
perturbations on the Oppenheimer-Snyder solution in this study. Thus we introduced, a dimensionless quantity 
related to the amplitude of gravitational perturbations, e ~ |<53/ii/|, an another one for the strength of the magnetic 
fleld, 77 ~ \B/{GAI^R^'^y^^\. Here Sgf^i, stands for the metric perturbation and B for the magnetic fleld strength 
while M and R denote the mass and radius of the star. Note that in this work we assume that the Lagrangian 
displacement of the fluid satisfles the condition ^ e. It should be emphasized that a perturbative treatment 

of the magnetic flelds is good enough to describe the strong magnetic flelds met in magnetars because dynamics 
of the stellar collapse is basically governed by gravity. In other words, the ratio of the magnetic energy £m to 
the gravitational energy £g is sufflciently small even for magnetars, i.e. a typical value of 8m I for magnetars is 
approximated by EmI^G ~ B^ /{GM^R^^) « 10^^ x (i3/10^^[G]). As for gravitational perturbations, this paper 
focus on the axial parity perturbations as the flrst step. The axial parity gravitational waves are treated with the 
covariant gauge-invariant formalism. We further assume that the two expansion parameters, e and 77, satisfy e ~ 77^. 
In other words, we consider the gravitational radiation directly driven by the magnetic fleld of the dust sphere. 

This paper is organized as follows. In ^|ll]we make a short introduction on the orders of the various perturbative 
quantities that we use, then we briefly describe the gauge-invariant perturbation theory, and flnally we show the form 
of the basic equations describing the magnetic flelds. Next, in jjllll we describe the formulation for the magnetized 
homogeneous dust collapse including an analytic description of the background spacetime, the magnetic flelds, the 
axial parity perturbation equations and the junction conditions at stellar surface. In i jlVl we describe the details 
of numerical procedures employed in this study, while the code test are presented in S|Vl The results related to the 
efflciency of the collapsing magnetized homogeneous dust spheres in gravitational wave emission are shown in ijVIl 
The final section. Will is devoted to discussion and conclusions. In this paper, we adopt the unit of c = G = 1, where 
c and G denote the speed of light and the gravitational constant, respectively, and the metric signature is (— , +, 
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II. PERTURBATION THEORY PRIMER 
A. Ordering of Perturbations 

Since the energy of the magnetic field is much smaller than that of the gravity even for magnetars, as mentioned in 
the previous section, it is reasonable to treat the electromagnetic field in the dust sphere as small perturbations. The 
perturbed metric, g^j/, the perturbed four- velocity of the fluid, u^, and the perturbed electromagnetic tensor, F^^, 
can be expanded as: 

gt,u = gi,u + Sgf,u + 0{e'^), (1) 
= ui' + Su'' + 0{e^), (2) 

= ^M'^ + -^i^M- + > (3) 

where 17^,^ is the background metric tensor, u^, the four- velocity of the fluid and F^^ the electromagnetic tensor. 
Both f^^jy and are defined as solutions of a collapsing spherical dust sphere in the absence of any electromagnetic 
field. For convenience we introduced two small dimensionless parameters related to strength of the magnetic field 
and to amplitude of the gravitational waves, i.e., rj ~ |i?/(GM^i?~'*)^/^ | and e ~ |<55piy|. Moreover, we assume that 
the Lagrangian displacement is small i.e. ~ e. Finally, we need to mention that in this study, we consider an 

infinitely conductive fluid , i.e. we make use of the so called ideal magnetohydrodynamic approximation. Thus, the 
master equations for describing the magnetic field are given by the perfect conductivity condition, F^nV^ = 0, and the 
Maxwell equation F^i^^q, -|- F^a./i + F^fi^u = 0. The fisrst order (ry^) form of the above two conditions will be written 
as: 

SF^.u" = , (4) 
(JF^^.c + (JF^cM + 5Fat.,v = , (5) 

which determine the magnetic field corrections to the spherical dust sphere. Up to this order of approximation, the 
variations induced by the presence of the magnetic field do not affect the spherical symmetry of the system, since the 
Lorentz force which induces deformations in the geometry is of second order (77^). 

In a similar fashion, both the Einstein tensor and the energy-momentum tensor can be expanded in powers of e 
and 77 as 

G^, = G^, + 6G^, + 0{e^), (6) 

r(r^-T(f) + 5r(f) + o(e2), (7) 

fl^.''^^5TlS''^+0{e^^), (8) 

where xj^j^^ and T^^^^^ stand for the energy-momentum tensors for the fluid and for the electromagnetic fields, 

respectively, while STj^^fj^^^ is of second order in 77. The Einstein equations of order rj'^e^ are the evolution equations 
describing the unperturbed spherical dust collapse. Here we focus in the study of the influence of magnetic field 
on the efficiency of gravitational wave emission during the collapse, and thus we consider only those terms of the 
approximation that will significant in this study. That is we omit terms such as rj'^e^ and we further assume that 
e ^ . In this order of approximation, the Einstein equations of order e are reduced to the following form 

5G^, = 8^{<5r^f ) + <5r^f *^)} + 0(e2) = S^5T^, + 0{e^) , (9) 

which describes gravitational perturbations driven both by the magnetic field and the fluid motions of the collapsing 
dust sphere. 

B. Gauge-Invariant Perturbation Theory 

The gauge-invariant perturbation theory for spherically symmetric background spacetime has been formulate d by 
Gerlach and Sengupta |23] while its covariant formulations has been developed by Gundlach and Martfn-Garcia [29l |. 
Here we only briefly describe this formalism for special case of axial parity perturbations. 
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1. Background Spacetime 

A spherically symmetric four dimensional spacetime Ai can be decomposed as a product of the form M = xS^, 
where A4^ is a 2-dimensional (1+1) reduced spacetime and a 2-dimensional spheres. In other words, the metric 
gfj^^ and the stress-energy tensor T^^ on A4 can be written in the form 

= diag{gAB,R^lab), (10) 

T^,=diSig{TAB,QR^lab), (11) 

where qab is an arbitrary (1 + 1) Lorentzian metric on R a. scalar on A^^, Q some function on and 7ab is the 
unit curvature metric on S^. Note that if the background spacetime is spherically symmetric then Q = . Here 

and henceforth the Greek indices denote the spacetime components, the capital Latin indices the components, 
and the small Latin indices are used to denote the S'^ components. Furthermore, the covariant derivatives on M, 
M'^, and are represented by ia, and -a, respectively. Finally, the totally antisymmetric covariant unit tensor on 
is denoted as Eab and on as Eat- 



2. Nonradial Perturbations 



As mentioned before, in this paper, we only consider axisymmetric axial parity perturbations both for the metric 
Sg^ii, and the matter perturbations 5T^^, which are given by 



5T^^ = 











(12) 
(13) 



where = e^^Pi-b while Pi stands for the Legendre polynomial. The gauge-invariant variables of the perturbations 
are then defined as 



kA = ht^' - h\A + 2hvA, 

La = AtT^^ - QhT'"\ 
L = At- Qh, 



(14) 
(15) 
(16) 



where va = R\a/R (29| . In terms of the gauge-invariant variables, the master equations for the axial parity pertur- 
bations are given by 



+ {l-l)il + 2)k 



^ - WttR^L^, 



(17) 
(18) 



{R'L^),A^{l-m + 2)L. 



(19) 



C. Basic equations for the magnetic field 

As mentioned earlier, the electromagnetic field perturbations, Ji^^jy, are governed by the Maxwell equations, i.e. 

SFf,,,^ + SF,,,f, + 6F,^,, = 0, (20) 
SFf"^^ ^ AttSJ'', (21) 



where SJ'^ is the perturbations of the current four-vector. Note that equations ([20|) and ([21]) are correct up to order of 
jy^e". The perturbation of the electromagnetic field energy-momentum tensor, dTlul^^\ in this order of approximation 
has the form 

1 / „ 1 , „ \ 

(22) 



K^""^ = ^ {^SP^^o.SF^pg''^ - ^g^JF^pSFx^g^^g'"''^ 
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The electric and the magnetic field associated with the four-velocity of the fiuid u"^ are defined as 

= dF^^u", (23) 
= ^e^^c.ffu'^SF"^ . (24) 

Finally, we remind that in this paper we consider infinitely conductive fluids, i.e. the ideal magnetohydrodynamic 
approximation has been adopted, according to which = SF^^u" = 0, where u" is the unperturbed four-velocity of 
the infinitely conductive fiuid. 

III. MAGNETIZED HOMOGENEOUS DUST COLLAPSE: FORMULATION 

A. Background spacetime for perturbations 

Here we briefly describe the background spacetime which will be later endowed with a magnetic fleld. We consider 
perturbations aromid a homogeneous spherically symmetric dust collapse described by the Oppenheimer-Snyder (OS) 
solution, whose line element inside the dust sphere is given by 

= -rfr^ + i?2(r)[dx^ + sin^ x{de^ + sin^ ed4)'^)] , (25) 
= B?{ri)[-d7f + dx^ + sin^ x{de^ + sin^ 0d<t>'^)] , (26) 

where % is a radial coordinate defined in the range of < x ^ Xo- Here Xi) is the stellar surface and it is assumed 
that xo < 7''/2. In the line element defined earlier, R{vi) is the scale factor and t{vi) is the proper time of an observer 
comoving with the fluid, defined in terms of the conformal time r] as follows 

M 

R{v) = -^{l+cosv), (27) 
sm'" Xo 

M 

t{v)^ —^iv + smr]), (28) 

sm Xo 

where M is the total gravitational mass of the dust sphere. The energy-momentum tensor for the dust fluid is written 
as 

T^f ) = pu^u,, (29) 

where p is the rest mass density given by 

o • 6 

and denotes the four- velocity of the dust, descibed in terms of comoving coordinates as 

u'^ = (5'; or u^ = R{ri)5^'^ (31) 

where 5'^^ means the Kronecker delta. The spacetime outside the dust sphere is described by the Schwarzschild metric, 
i.e., 

ds'^ = -f{r)de + f{r)-^dr'' + r^{de^ + sin^ Odcj)''), (32) 

where /(r) = 1 — 2M/r. From the junction conditions at the surface of the dust sphere, we obtain the relationships 
between the {rj, x)-coordinates and the {t, r)-coordinates, given by 

= i?(??)sinxo, (33) 



[(r,o/2M)-l]i/2 + tan(r//2) 



[(r,o/2M)-l]i/2_tan(77/2) 



where Vso = rs{t = 0) = 2M/ sin^ xo is the initial stellar radius in Schwarzschild coordinates. 
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B. The Magnetic field of tlie star 



As mentioned earlier, we consider weakly magnetized dust spheres in which the magnetic effects on the dust fluid 
are treated as small perturbations on the OS solution. Moreover, for the sake of simplicity we assume that the 
electromagnetic fields are axisymmetric. Thus, perturbations of the electromagnetic fields, SF^^^, and the current 
four-vector, ^ J^, can be described in terms of the Legendre polynomial by the following relations 

(35) 
(36) 
(37) 
(38) 
(39) 
(40) 

(("'S-i") . (41) 



SFo3 = 


-SF30 


= ei sin 9 dg Pi 


5Fi3 = 


-6F31 


= bi sinOdePi 


6F23 = 


-6F32 


= 62 sin 9Pi^j . 


SFoi = 


-SF,o 


= e2-P/„ , 


SF02 = 


-SF20 


= esdgPij^j , 


5Fi2 = 


-6F21 


= bsdePi^r - 


5J^ = 





Notice, that here we have used Im to denote the angular quantum number with respect to the electromagnetic fields 
to discriminate it from the one for the gravitational waves I. 

In the interior of the dust sphere, the perfect conductivity condition SF^i^u'^ = is reduced into SFq^ = 0. This 
assumption leads to the following simplifications 

61=62-63 = 0. (42) 



By direct substitution of equations ([35]) through (j40|l into the Maxwell equation ([20|) . we obtain the basic equations 
describing the magnetic fields, which have the following simple form 

a^foi = d^b2 = d^bs = 0, (43) 
lMilM + l)bi+d^b2^0. (44) 

The first of theses relations, equation (j43p . suggests that a comoving observer does not observe any change in the 
magnetic field distributions. All the components of the electromagnetic fields can be determined through equations 
(j42ll . (|43|) . and (|44p . The Maxwell equation ([2T|) that we still have not use can be regarded as the definition of the 
current four- velocity. This implies that the perturbations of the current four- velocity can be written as follows: 

h = 0, (45) 



Im{ Im + 1) 

47r i?2 sin^ X ' 



.'" = -^, (47) 

The electromagnetic fields outside the star are also given by similar expressions to equations (l35 |) -(|40 |) . In order to 
avoid mixing of the various quantities inside and outside the star we will indicate the ones in the exterior with a tilde i.e. 
ei, 61, and so on. Since the exterior of the star is vacuum, we cannot make use of the perfect conductivity condition 
there. Instead, we make an alternative assumption that is we demand the vanishing of the current perturbations 
outside the star i.e. = 0. This assumption simplifies considerably the Maxwell equation ((2T|) leading into thre 
following set of equations 

drSr^e2)-lM{lM + l)i3^0, (49) 

dt{r^S2)-lM{lM + l)fb3^0, (50) 

dtii-dr,{fbi)-^b2^0, (51) 

where is the tortoise coordinate defined as = r -I- 2Mln(r/2M — 1). Furthermore, from equation (|^n|) . we get a 
set of equations similar to equations and (M)) . given by 

62 - dri3 + dth = 0, (52) 
lM{lM + l)ei+dtb2^Q, (53) 
/m(W + 1)61 4-9^62 = 0. (54) 
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The six equations for the electromagnetic fields, i.e. equations p9|) through (|54p. can be reduced to two decoupled 
wave equations 



-dTh^ + dlh - ^J^i^^fh, = , (55) 
-dKr^-e^) + dl (r'e^) - ^^^^^^/(r^e-^) = . (56) 



Moreover, these two wave equations can be rewritten in terms of the double null coordinates, u = t — r^ and v = t + r^, 
as 

dudv 4r^ ' 

At the surface of the star, we implement the following junction conditions for the electromagnetic field 

n^B^ = n^B^, (59) 
= q^^E,, (60) 

where n^, n'' are the unit outward normal vector to the stellar surface defined in the interior and the exterior 
coordinates, respectively, while , q^^ are the corresponding projection tensors associated with and fi'^. Therefore 
the junction conditions reduced to the following set of relations 

b2 = b2, ei + —bi =0, 63 + -77763 = 0, (61) 

where 

— = ^/tanxo- 62 

Concluding we mention that in this paper we focus only on dipole electromagnetic fields, i.e., electromagnetic fields 
associated with Im = 1- Observations actually are in favor of the existence of dipole electromagnetic fields and 
moreover these fields can drive the quadrupole gravitational radiation as we will see in the next section. Finally, the 
formalism developed here accounts for electromagnetic fields which lie both inside and outside the star, in this study 
as a first step, we take into account only magnetic fields confined in the stellar interior. We therefore assume that 
Ci — bi — foi i — 1, 2, 3. 



C. Basic equations for the axial parity perturbations 

1. Interior region of the star 

As we mentioned earlier we will use the gauge-invariant formulation in the treatment of perturbations of the OS 
spacetime. The gauge-invariant form of axial perturbation equations (|17p and (jlSp for the OS spacetime is reduced 
to the following set of equations 

+ d^k^ - WnR^L = 0, (63) 
d^iR^Usin'^ x) + il-l)il + 2)^- 167rL,, sin^ x = 0, (64) 

^^^^(i^'nsin^ x) + {l-l){l + 2)^- IGttL^ sin^ x - 0, (65) 

where 11 is the gauge-invariant quantity, defined as 

1 



n - , 

i?2 



i?2sin^X/ V^^sin^X 



(66) 
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The regularity condition at the stellar center, suggests the introduction of a new function 11 defined as 

n = (i?sinx)'-'n, (67) 

which is analytic at the stellar center. By using equations ([64|) and ([65]) . one can derive a single wave equation for n, 
given by 

where -R(O) — 2M/ sin'^ xo- Moreover the equation of motion p9p is rewritten as 

- dniR^Lr, sin^ x) + d^{R^L^ sin^ x) ^ {1 - + 2)R^L. (69) 
The perturbation of the energy-momentum tensor ST^i, can be splitted into two parts as shown before: 

5T,,. = ^r(f) + JT(f*^), (70) 

where T^t^^ and T^^*^'* are the energy-momentum tensors for the dust and the electromagnetic field, respectively. 
Since axial parity perturbations of the four- velocity of the fluid, defined as 

Su^ = [0,0, P{T,x)Si), (71) 

the expansion coefficients of StI;^^^ introduced in equation are given by 

Atf)=/3pz.„ (72) 
Ai(*^) = At(*^) = 0. (73) 

In this paper we constrain our study to the quadrupole gravitational radiation emitted by axial parity perturbations. 
The reason is that quadrupole radiation directly couples with dipole magnetic fields and it is dominant component 
for gravitational wave emission. As discussed before, the perturbations of the energy-momentum tensor for the 
electromagnetic field, 5tIv^'^\ are of order ~ ((5F^j/)^. Therefore, we cannot achieve separation of variables for 
<5r^^*^' even if the background spacetime is spherically symmetric. Since we assume dipole electromagnetic fields 
{Im — 1), then ST^^^^^ will contain terms associated with I — and / = 2. For the detailed calculations of various 

components of STjlfl^^^ inside the star, see Appendix 1X1 Finally, the expansion coefficients for Sxjlu^^^ associated with 
I — 2 are given by 



Atf = 0, (74) 



At(^^^) = lUl- (76) 



Thus we can derive the gauge-invariant quantities for the total matter perturbations (the dust fluid and the magnetic 
field). La and L, which have the following form 



L,^ = -Rpp, (77) 

127ri?2 sin^ X 
6163 



Lx^^TWZTW^.^ (78) 



127ri?2 ■ 

Substituting equations (|77p through ((79)) into equation ([BO]). we get the following equation of motion for /3 



127ri?3psin^X 



(79) 



(80) 
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notice that in this derivation we have used the perturbed Maxwell equation (|44)) . Using the relations, R^p = 
3R{0)/8tt — const., drp2 — 0, and 9^,63 = 0, we can analytically integrate equation ([50)1 with respect to confor- 
mal time 77, to get the solution 



2V 



9i?(0)i?'+i sin'+-* X 



(/-!)(/ + 2) 



i+i 



where (3 = sinx)'+^ and /3o(x) is the initial distribution of /3. Following Ref. [2E 
different definitions for the initial distribution /3o(x) i-^- 



/3o(x) = lli{= const.), 

^i?(0) sinx 



/3o(x) =Z^3Cxp 



Rc 

/ i?(0) sinx - ^^0 
V Rc 



/3o(x), (81) 
in this paper, we adopt three 

(82) 
(83) 

(84) 



where Ui is an arbitrary constant and Rc is a scale factor describing the inhomogeneity of the fluid velocity's initial 
deistribution. Here following Ref. [1^ we choose Rc = rso/3. 

In the actual numerical calculations, we use the two null coordinates, {u = rj — x and v = rj + x) and the master 
equation (j68p is rewritten in these coordinates as 

n = 5(77,x), (85) 
(86) 



d^n + l /^cosx , dnR\ dll l + l fcosx driR\ dll ^ {21 - + 2)R{0) 



dudv 



smx 



R du 



smx 



R dv 



S{r,, x) = 47ri?V {(^ + 1) cosx/3 + smx{dj)} 



8R 

2b2b3{d^R) 



3i?'+3 sin'+2 



X 



In summary, inside the star, our basic equations describing the gravitational perturbations are equations (I85p and 
(I86p . The source term S{r],x) is given uniquely by the function P{rj,x)i shown in equation (j81[) . for a given initial 
distribution of the electromagnetic field and the fluid velocity perturbations. 



2. Exterior region of the star 



The Oppenheimer-Snyder solution (|26|) for the interior is matched with the Schwarzschild solution (|32|) for the 
exterior and the master equations for perturbations (fTT)) and in the interior reduce to the well known Regge- 
Wheeler equation in the exterior, which has the form 



V{r) 



where $ is the Regge- Wheeler function, related to the gauge-invariant variable 11 through the relationship 




(87) 

(88) 

(89) 



If the electromagnetic fields do not vanish outside the star, we need take into account their influence on the gravita- 
tional radiation cmmitted during the collapse. In the exterior the only non- vanishing contribution to the perturbations 
of the energy-momentum tensor ST^^i, is the one from the electromagnetic field. The expansion coefficients for (^T^^ 
associated with the I = 2 axial parity perturbations are then given by the following formulae 



1 

12^ 
1 



12^ V/ 



fe2bi 
1 . . 



-eiC2 



-46263 



1 



127r V/ 



1 



:eie3 - 



(90) 
(91) 
(92) 
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Actually, in Appendix |B1 we present the detailed form of the various components of the perturbed energy- momentum 
tensor for the electromagnetic field outside the star. 

By using the vacuum Maxwell equations (|49|) -((54 |) . we can rewrite Lt, Lr, and L in terms of 62 and 62 as 



247r 



62(3^62) 6362 - &2(3re2) 

r 



62(9*62) - 62(9*62) 



L = 



487r 



24^/ 

2ri2{dtb2) + r\dtb2){dri2) - r\dt~e2){drb2) 



(93) 
(94) 
(95) 



Finally, the Regge- Wheeler equation ([57]) can be rewritten in terms of the double null coordinates, u = t — r^, and 

V — t + r^,, as 



S + Inr)^ = ^ [r'S.idrM - 62(9.. r2e2) 



(96) 



where we have modified the right-hand side of the Regge- Wheeler equation by using equations (|55|) and ((56|l . Still, 
since in this study we don't take into account the influence of the electromagnetic field outside the star, the right-hand 
side of equation ([M)) vanishes. 



D. Junction conditions at the stellar surface 



In order to ensure that the spacetime is regular at the stellar surface (x — Xo)i we impose three junc- 
tion conditions for the case of axial parity perturbations; first we demand the continuity of 11, second that 
n"^n|^ — 167r(i?sinx)~^M"^-^A = n^^\A — l&nr~^u"^LA, and the last is u"^n|^ — u^YI^a- These boundary con- 
ditions arise from the continuity conditions for the induced metric and the exrinsic curvature [29j . Therefore the 
junction conditions are explicitly given by 



n = n, (97) 

-Z + W + ^^^2^ = ((9,,i?)sinxo - -Rcosxo) ^ + ((9,,i?) sinxo + ^cosxo) 

Rsm Xo J J 

--^ (^LtRcosxa + fLr{driR)smxoj , (98) 

Z W 

Z + W = {Rcosxo - {driR) sin Xo) y + (-R cos xo + {dr,R) sin Xo) y , (99) 

where Z = dH/du, W = dW/dv, Z dn/dii, and W = dil/di. The last two conditions, jMl) and jSSl), can be 
rewritten as 

W = (i?cosxo + (9^^)sinxo) - 4^ (Zti^cosxo + /ir(9,,i?) sinxo) - ^^^2 ' (1^0) 

/ /? ^ ^ -Rsm^Xo 

. Z S-K { ~ ~ \ \ Stt/Sp 

Z = (i?cosxo - (9^^)sinxo) T + (Lti? cos xo + /Lr(d,,i?) sinxo + „ . 2 — ■ (101) 

J jr-' \ / i?sm Xo 



IV. NUMERICAL PROCEDURE 



In this section we describe the numerical procedures that we will follow and the way that we generate initial data. In 
order to simplify the numerical procedure and to set the initial data both in the interior and exterior of the collapsing 
configuration, we divide the background spacetime into three regions named I, II, and III, as illustrated in Fig. [1] 
Region I represents the stellar interior, while regions II and III the exterior spacetime. Region III is separated from 
region II via the null hypersurface defined by v = vq, which is generated by the ingoing null rays emitted from the 
point where the stellar surface reaches the event horizon. Note that it is sufficient to consider the regions I, II, and 
III because the gray area in Fig. [1] is causally disconnected from the stellar interior at 77 = when the magnetic 
fields confined inside the star. To solve the wave equation numerically, we make use of the finite differencing scheme 
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proposed by Hamade and Stewart [stI, in which the double null coordinates (u, v) are employed in region I and {u, v) 
in regions II and III, respectively. Notice that we integrate the wave equation in region I by using a first order finite 
differencing scheme to avoid numerical instabilities appearing near the stellar center, while Hamade and Stewart's 
original scheme is of a second order finite differencing scheme. 




f] = >> 
stellar surface 



FIG. 1: A schematic description of the Oppenheimer-Snyder spacetime for the collapsing model in characteristic coordinates. 
Region I denotes the stellar interior while regions II and III correspond to the exterior. The stellar surface, where r = or 
X = xo, is the boundary between regions I and II, and the stationary region outside star is indicated by gray. 



A. Initial Data 



In order to initiate the numerical calculations, we need to provide a data set on the initial hypersurface for the 
quantities fl, du^, dyll, (3o, 62, and 63 for the stellar interior, and $, du^, and dy^ for the stellar exterior. Following 
[22l |. we assume that the initial perturbations are "momentarily static" . Outside the star, the momentarily static initial 
condition for the metric perturbations is given as a static vacuum solution of equation (j87|) in terms of hypergeometric 
functions 

- ^(7^ (^) ' - M + 3, 2/ + 2; ^) , (102) 

where qi is a constant representing the multipole moment of the star. Here we assume that qi — 2M. We remind that 
there is no electromagnetic field outside the star. Since it is a static solution, the initial perturbation outside the star 
(|102p does not evolve until a light signal from the stellar interior arrives there. Finally, set the initial data i.e. the 
static solution (|102p . on the null hypersurface u = uq for the characteristic initial value problem. 

As for the initial condition inside the star, we have to give a data set not only for metric perturbations but also for 
fluid and electromagnetic perturbations. As mentioned earlier, the initial fluid distribution f3o is given by equations 
([82)1 - (|84|) . while the initial distribution for the magnetic fields 62, &3 will be discussed later in WII Finally, the 
momentarily static initial data for metric perturbations ft in region I defined by the conditions 9,^11 = and d^fi = 0. 
The data sets at 77 = for du^l and dyll are then defined via the relations 9„II = — (9^n)/2 and dyll — {dy^Il)/2, 
respectively. The momentarily static distribution of 11(77 = 0) due to the conditions 9^11 — and d^Il = 0, is given 
as a regular solution of 



The regular solutions of equation (|103p must be smoothly connected to the static exterior solution (|102p through the 
stellar surface. This leads to a boundary condition for II(?7 = 0) at the surface of the star described by the following 
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equation 



21 + 1- ^\U + tanxon,, + — + . * , , = 0, (104) 



where Fi is the abbreviation for the hypergeometric function Fi{l — 1,1 + 3, 21 + 2; 2M/r). FinaUy, the regularity at 
the stellar center requires that the function 11 is analytic for x ^ 0. 



B. Boundary Conditions 



The boundary conditions for the numerical integration are the regularity condition at the stellar center and that 
there are no incoming waves at the infinity. The regularity condition at the stellar center demans that d^Ii = 0, which 
is reduced to 

ou ov 

Finally, for the no incoming radiation condition at the infinity, we adopt the condition d^/du = (see, e.g., [37|V 



C. Special Treatment of the Junction Conditions near the Event Horizon 

When the stellar surface reaches the event horizon, the junction conditions discussed ealier in i jlll Dl cannot be used 
any more because the terms related to /^^ diverge. Instead of these junction conditions, following 2^, we impose 
the following junction conditions on the null surface oi v = vq in the vicinity of the point Ti. in Fig. [TJ 

ttEH _ ttJV,„,, 

n-n^- + {r-r---), (106) 

n"" - n-»^^ + - r---) , (107) 

where IT" and r" are the values of 11 and r on w = vq at n-th time steps, while A'max denotes the total number of 
time steps in region II, and r^^ — 2M. 



V. CODE TESTS 



In order to check our numerical code, we have calculated the quadrupole gravitational radiation emitted during the 
collapse of a non-magnetized homogeneous dust sphere (perturbations of the Oppenheimer-Snydcr solution), which 
has been already studied by several authors, e.g., [Ij,!!^. For the test we will consider the collapse of the homogeneous 
dust sphere which is initially at rest. Therefore, we have to provide the initial radius of the dust sphere to begin the 
numerical integration. Since the amount of gravitational radiation emitted during the collapse of a non-magnetized 
homogeneous dust sphere is the "typical value" with which we will compare the energy emitted during the collapse 
of a magnetized homogeneous dust sphere, we will briefly summarize these results. 

In the present calculations, the number of the spatial grid points inside the star (region I) is chosen to be = 1000. 
Using this number of grid points we manage to obtain numerical solutions and results with acceptable accuracy. In 
region III, the step-size for integration is determined by the relation Au — (wmax — uq)/Nu, where Wmax = imax — r^oh- 
Here, imax is the expected maximum time for observation and r<,ob is the position of the observer described in tortoise 
coordinate units. In this paper, we assume that i,nax = 2000M while the fiducial observer is at rob = 40M. 

As a first step we confirm the convergence of our numerical code. For this purpose, by varying the value of N^, we 
calculate the total energy radiated in gravitational waves during the collapse which will be characterized by the initial 
radius of the dust sphere, r^o = 8M, and the /3o = const, velocity distribution. The radiated energy is estimated 
by integrating the luminosity of gravitational waves i-cw,; with respect to time. Here, the luminosity Lgw,i of 
gravitational waves is defined by the relation (see, e.g., [22, 25]) 
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The outcome of this calculation i.e. the energy emitted in gravitational waves during the collapse as a function of Nu 
is shown in Fig. [5] and tabulated in Table HI From Fig. O we conclude that the amount of the total energy emitted 
during the collapse converges for > 10000, thus we assumed in all numerical runs = 10000. In Fig. [3l we give 
the relative error in the total emitted energy obtained by our numerical code as a function of At2. Here, the quantity 
— {E{Au) — Em)/E„i stands for the relative error, then E{Au) denotes the total energy emitted for various values 
of Ail and Em is the energy for some maximum value Nu — 20000. It is obvious from Fig. [3] that our numerical code 
achieves second order accuracy. 



6 X 10-' 



a 

u 



2 X 10-^ 1 ^ ^ — ^ 1 

300 1000 10000 30000 

FIG. 2; The total radiated energie of quadrupole gravitational waves from the OS collapse characterized by the initial radius 
of the dust sphere r^o = 8M and /3o = const, velocity distribution, as function of the number of grid points in region III. 



TABLE I; The total radiated energies of quadrupole gravitational waves from the collapse characterized by the initial radius 
of the dust sphere rao = 8M and the /3o ~ const, velocity distribution. 



Nu Au emitted energy / (2M) 

500 3.929M 5.614 x 10"^ 

700 2.806M 3.641 x 10"^ 

1000 1.964M 2.964 x 10"^ 

5000 0.3929M 2.255 x 10~^ 

10000 0.1964M 2.197 x 10"^ 

20000 0.09822M 2.175 x 10"^ 



Next, let us compare the total energy emitted during the collapse as it has been calculated by our numerical code 
with the results by Cunningham, Price, and Moncrief 22] (CPM1978). In Fig. [H we show the total energy emitted 
in gravitational waves during the collapse as function of the initial radius of the dust sphere r^o- In this figure, the 
results of CPM1978 are indicated by the fiUed-squares, while the other symbols represent our results. The results for 
the initial distribution of the fluid velocity defined by equation are indicated by circles, those defined by equation 
([55)1 are indicated by triangles while for the distribution defined by equation ([5^ we used squares. Finally, the results 
obtained by using a numerical code with the first order accuracy, are indicated by the gray asterisks. This figure 
shows that there are small differences between our results and those obtained by CPM1978. We however observe 
that the results of CPM1978 agree well with those with the first order accuracy (compare the fiUed-squares with the 
asterisks). Therefore, we conclude that our results are in quite good agreement with the results of CPM1978. 
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r,o = 8.0M 
Po = const. 
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FIG. 3: Convergence test of the numerical code. The vertical axis denotes the "relative error", that is the ratio of E{Au) — Em 
over Em, where E{Au) is the emitted energie for various of Au and Em is the energy for A''^ — 20000. The dashed line is 
cx (Ait)^ suggesting that as Au becomes smaller the "relative error" reduces as (Au)^. 
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A A ot exp[- (i?(0) sin x/R,f\ 
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FIG. 4: The energy emitted in gravitational waves from the homogeneous dust collapse without magnetic field as a function of 
the initial stellar radius where r^o ~ 8A4, 12M, 16M, and 20M. The fiUed-squares correspond to the results by 22] while the 
rest correspond to our results. The different marks corresponds to different initial distribution /3o ■ The asterisks correspond to 
numerical results taken by a first order code. 
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In Fig. \E\ we show the waveforms of the quadrupole gravitational radiation from the collapse with an initial radius 
Vso = 8M and for three different initial distributions of the fluid velocity defined by equations ([5^ . and ([M)) . 
In this figure, the left panel displays the waveforms as functions of the time, while the right panel displays the absolute 
values of the amplitudes as functions of the time in a log-log plot. Fig. [5] shows that the first part of the waveform 
is characterized by the quasinormal ringing while at the late times follows a power-low tail, as found in [2^ . (For a 
review of quasinormal modes for compact objects, see, e.g. (ssj.'l From the waveforms, we estimate the frequency 
of the fundamental quasinormal mode to be 2Muj — 0.746 -I- 0.179i, which agrees very well with the quasinormal 
mode frequency estimated by Chandrasekhar and Detweiler [s^. As for the late-time tail of the gravitational waves, 
in the right panel of Fig. O we find that the amplitude decays as t~'^ at late times, which is in good agreement 
with the analytical estimate of Price \M, that is i-(2'+3). The accuracy in the estimates of the quasinormal mode 
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frequency and the late-time tail therefore suggest that our numerical code is accurate and reproduces all previously 
known results. 

Besides the tests of the code the following basic properties of the gravitational radiation emitted during the collapse 
of the non-magnetized homogeneous dust should be emphasized. First, as it shown in Fig. [H the total radiated 
energy does not critically depend on the distribution of f3o, this has been also observed in CPM1978 and second, 
as shown in Fig. [5l that small modulations appear just after the onset of the collapse only for the case of f3o oc 
exp[-(i?(0) sinx/i?c)2]. 




t [M] t [M] 

FIG. 5: The waveform of the quadrupole gravitational radiation emitted during the collapse of the non-magnetized homogeneous 
dust, as function of time. The initial radius of the dust sphere is set to rso = 8M while the fiducial observer is set at r = 40M. 
In the right panel the amplitudes of the gravitational waves are shown in a log-log plot and the late-time is compared with 
its theoretical value Three initial distributions of the fluid velocity f3o were adopted and one can hardly trace their 

influence on the waveform. 



VI. GRAVITATIONAL RADIATION FROM THE COLLAPSE OF THE HOMOGENEOUS 

MAGNETIZED DUST SPHERE 

A. Initial distribution of the magnetic field and magnetic eflFects on the gravitational radiation 

For the calculation of the gravitational waves emitted during the collapse of a magnetized dust sphere, one needs to 
provide the initial distribution (profile) of the magnetic fields, i.e., to set up the functional forms of 62 and 63 on the 
hypersurface 77 = 0. In practice, one has the freedom in choosing the initial distribution of 62 and 63 as the following 
two conditions are satisfied: (i) the regularity condition at the stellar center, (ii) the junction condition (pTj) at the 
stellar surface. Since here we made the assumption that the magnetic field is confined inside the star, the second 
condition is reduced to 62 (xo) = 0. In this work we assume a dipole magnetic field, and for this specific geometry it 
is convenient to introduce two new quantities for the description of the magnetic field, &2 and 63, defined as 

hix) = iRsinx)^b2{r],x), (109) 

b3ix) = iRsinx)%iv,x)- (110) 

If the new variables 62('7,x) and b^{r],x) are analytic at x = 0, the regularity condition at the stellar center for the 
magnetic field is automatically satisfied. Then for quadrupole perturbations (/ = 2), the source term in the wave 
equation, S{ri,x), and the perturbation for the four- velocity Pirj^x) are given by the following expressions: 

S {r^, X) ^ R^P{W cos x+{d J) sinx] - ^^^^^'"""^ (IM^ B^B^hohosinx, (111) 
Piv,x) = B2B3 {620630 tanx + 620(9^630) - {dMho} + (^^) Mx), (112) 
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where &20: ^30 are dimensionless functions of Xi defined as 

^2(?7,x) 
h{Tl,X) 



m 



(113) 
(114) 



Here B2 and B3 are arbitrary constants related to the strength of the magnetic field. It should be emphasized that, 
as shown in equations (|llip and (|112p . all terms related to the magnetic fields in the source term S{ri, x) vanish when 
ry = 0. In other words, the magnetic field does not affect the momentarily static initial data for metric perturbations. 
The geometry of the magnetic fields when the collapse sets in is practically unknown. Based on this freedom we 
adopted the following two types for the initial distribution of the magnetic field: 



(I) : 62o(x) = l-2(^ 
(II) : &2o(x)=^3o(x)=4 





2L. 




Xo 






i- 


r- 


Vxo 





63o(x)-l- - 

Vxo 

4" 



(115) 
(116) 



Note that the maximum value of &20 and 630 are chosen to be one in the interval < x < Xo- For the first profile 
function (I) the magnetic field is stronger in the center of the sphere while for profile function (II) the field becomes 
stronger in the outer region. 

It is more convenient to explore the effects of the magnetic field on the efficiency of the gravitational radiation 
emission during the collapse, if one introduces the dimensionless parameter a, defined as 



Then the source term in the wave equation, S'(??,x)j can be split as follows 



where 
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S(B) 



47ri?(0)V 
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Stt 



9i?(0)2 
sin 77 



{3/3oCOSx+ (9x/3o)sinx} , 



sm X 
620^30 sinx 



dx\ ^20^30 tan X + &2o(9^&3o) - (9x&2o)&30 sin x 



(117) 



(118) 



(119) 



(120) 



Note that S^f^'^ can be attributed to the incompressible fluid flow while S"^^^ to the Maxwell stress in the magnetized 
dust sphere. The splitting of the source term introduced with equation (jllSp suggests that any solution of the wave 
equations (|85|) and ([96)1 can be expressed as a linear superposition of the two solutions, ^'^^^ and that are 

independent of a, i.e. 



(121) 



where is the solution in the absense of magnetic field (a — 0) and is the solution for the case of 5 = S*^^-* 
i.e. when the gravitational field is initially stationary. Since we assume that the initial profile of the gravitational 
perturbations does not dependent on the existence of magnetic fields, the initial values of were set to zero. 

These assumtions, i.e. the splitting introduced by equations pisp and (|12ip . suggest modifications in the form of 
equation (jlOSp describing the luminosity in gravitational waves which gets the following form: 



1 



GW,l 



{(ci.(f)^+a^(ci>:r)^+M$:r)('i'r)}- 



s(S)^2 



167r (l-\){l + 2) 

While the radiated energy in gravitational waves during the collapse is then given by following relation 



(122) 



(123) 
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where Eq^^ and i?^^; stand for the total radiated energies associated with the solutions l>'^''^ and respectively, 
and Ccwd is an integral quantity defined by the product of ^'^'^^ and Their detailed form is: 

= T^TT^^^TTT^ / i^frdu, (124) 
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Cg».j = ^ / (*!?)(*!?)<'"■ (126) 



where Wmaa; is the maximum value of v. 

It is worth mentioning the following issue emerging from the study of equation (jl23p . It can be easily proved 

that since Egw,i is a quadratic function of a and E^y^i > then by definition, Eqw,i gets its minimum value, 
^Gw? ~ ^GWi ~ ^GW,i'^ /E^^i, for a = ac, where ac = —Cgw,i/E^^i. In other words, phase cancellation between 
the two components and of the gravitational perturbations become maximal around a — ac- As discussed 
earlier in fJVl E^^^ does not highly depend on the initial profile of the fluid velocity. Thus, the total radiated energy 
Egw,i practically depends on a and the initial radius of the dust sphere, Tsq, but not on (3q. 



B. Numerical results for the initial magnetic field profile (I) 

The numerical study for the influence of the magnetic field on the gravitational wave output during the OS collapse, 
we have mentioned earlier, that we used two quite different initial profiles for the magnetic field. These two profiles 
described by the equations (|115p and (|116p and represent magnetic fields which have their maximum either at the 
stellar center or near the surface. 

We first consider the case (I), given by equation (jllSp . As discussed in the previous subsection, the fundamental 
quantities for estimating the gravitational radiation from the collapse of the magnetized homogeneous dust sphere 
are ^^^^ and It is obvious that we can obtain solutions of equations ([55]) and (IM)) for any value of a in terms 

of l>'^''^ and l>(^) through the relationship (|12ip . In Fig. [SI we show the waveforms and l>(^) of the quadrupole 
gravitational radiation from the collapse characterized by an initial radius = 8M and an initial distribution of 
the fluid velocity /3o = const. We observe in Fig. [6] that the two waveforms, l'(^) and are almost in phase. 

Therefore, we expect the following properties of the gravitational wave amplitude: (1) for a > 0, as a increases, the 
gravitational wave amplitude $ increases monotonically and the phase of $ does not change, (2) for a < 0, as a 
decreases, the amplitude and phase of the gravitational wave amplitude $ show a more involved behavior due to the 
phase cancellation between the two gravitational wave amplitudes and It has been also found that the 

amplitude of the emitted gravitational waves is almost independent from the functional form of (3^ in agreement with 
our previous results for the non-magnetized collapse (see Fig. [5]). 

In Figure 7 we show the waveforms for a number of positive and negative values of a i.e. for a = —9, —6, —3, 3, 
6, and 9. We also show, in every panel, the shape of the waveform in the absense of magnetic field (a — 0) with a 
bold line. This figure verifies the previosly refered theoretical estimations i.e. that for a > 0, the phase of the various 
waveforms is almost the same as that of a = and the amplitude becomes larger as a increases. The amplitude of 
the waveforms is considerably smaller for a < 0, due to the phase cancellation effects, and decreases for smaller values 
of |a| on the other hand the phase shift is significant for large values of \a\. Finally, one can easily observe that the 
dependence of the amplitude on the initial profile of /3o is very weak. Moreover, we observe another effect related to 
the magnetic field i.e. there is a wave packet the actual quasinormal ringing. The appearances of this wave packet 
can be only attributed to the functional form of 

As it has been argued earlier, the total energy emitted in gravitational waves during the collapse of a magnetized 
dust sphere for any value of a can be calculated as a proper combination of the quantities E^^^, Eq^^, and Cgw,i 

via equation p23p . The values of Eq^^^ ^(^2^ ^'^'^ Cgw,2 with f3o = const, for four different initial radius of the 
dust sphere (rso = 8M, 12M, 16M, 20M) are summarized in Table [Til By studying equation (|123p we observe that 
the total energy emitted in gravitational waves has a minimum for collapsing models with a = ac < 0. Since _E^(^2's 
are almost independent from the initial distribution of the fluid velocity /3o, as shown in Fig. 31 the total emitted 
energy hardly depends on f3Q. These features can be seen in Fig. [5] where the total energy emitted in gravitational 
waves from the collapse is studyed as function of a 
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FIG. 6: The amplitude of quadrupole gravitational waves emitted during the collapse. The continuous line stands for the 
function and the dashed for The efficiancy of the collapse depends on the initial radius of the dust sphere Vao = SM, 

the initial profile of the fluid velocity /3q = const., and the initial pofile of the magnetic field, here is the case (I). The waveforms 
are monitored at r = 40M. 
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FIG. 7: Waveforms for gravitational radiation emitted during the collapse of of a magnetized homogeneous dust sphere with 
initial radius Vso = 8M and initial distribution of the fluid velocity /3o = const.. The gravitational waveforms are monitored at 
r = 40M. The thick continous line corresponds to the waveform of the non-magnetized collapse (a = 0). 
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The most important conclusion that can be drawn from this figure is that the total energy emitted in gravitational 
waves from the magnetized dust collapse can be about eleven times higher than the energy of the non-magnetized 
collapse for a = 10 and about five times smaller for a = olc- Thus, the effect of the magnetic field in the gravitational 
wave outcome during the collapse can be significant and might improve the possibility of detecting gravitational waves 
from this type of sources. 



TABLE II: The values of Eq^j, E^^i, and Cgw,i for quadrupole gravitational waves {I = 2) emitted during the dust collapse 
for four values of the initial radius of the dust sphere (vso = 8M, 12M, 16M, and 20M) while we assumed that /3o — const. 
The value, ac = —CGw,2/E^^2y i-^- the minimum of the emitted energy Egw,2 is also shown. 





rso = 8M rsi) = 12Af r^o = 16A'/ r,o = 20M 


4«,2/(2M) 
4^^_2/(2M) 
Cgw,2/{2M) 
ac 


2.197 X 10"^ 8.319 x 10"** 4.354 x 10"* 2.676 x 10"® 
1.296 X 10"** 7.527 x lO"'' 5.557 x 10"^ 4.472 x lO"'^ 
4.722 X 10"* 2.366 x 10"* 1.501 x 10"* 1.066 x 10"* 
-3.64 -3.14 -2.70 -2.38 
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FIG. 8; The total energy emitted in gravitational waves {I — 2) during the homogeneous dust collapse with magnetic field as 
function of a. The horizontal axis denotes the value of the parameter a representing the strength of the magnetic field, while 
the initial radius of the sphere assumed to be r^o = 8Af. 



Next we examine the dependence of the total energy emitted during the collapse on the initial radius r^o of the dust 
sphere assuming again that the initial profile of the fluid velocity is /3o = const. In Fig. [9l the waveforms from the 
collapse with an initial radius Tso = 20M are shown. From this figure, can be easily seen that the effect of the 
magnetic field becomes more pronounced as the initial stellar radius increases while the basic properties (dependence 
on a and phase shift) are similar to those of the rso = 8M case. The influence of the magnetic field in the gravitational 
wave output increases with increasing radius, a natural explanation is that the longer the collapse lasts the longer 
the magnetic field will influence the dynamics of the collapsing dust. Actually, if the collapsing sphere has an initial 
raidus r^Q — 8M then it takes t ~ 50M until one observes the first peak of the quasinormal ringing, while it takes 
t ^ 140A'/ for the r^o — 20M model. Additionally, we can observe that a wave packet before the quasinormal ringing 
is suppressed, compare with Fig. [51 

In Figure 1101 we show the total energy emitted in gravitational waves as a function of a for several values of the 
initial radius of the dust sphere. From this figure, we can see that the critical value increases as rso increases, see 
also Table nil while the emitted energy decreases as Vso increases. However the effect of magnetic field is stronger for 
the model with larger initial radius, i.e., the emitted energy varies with the strength of the magnetic field. Actually, 
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as we mentioned earlier the ratio of the energy of the magnetized cohapse over the energy of the non-magnetized one 
varies from 0.22 (a = ac) to 11.2 (a = 10) for r^o — 8M. For r^p = 20Af the same ratio varies from 0.050 (at a = ac) 
to 25.7 (at a = 10). 




FIG. 9: Waveform of gravitational waves for I — 2 for (a) a > and (b) a < 0. The initial radius is Vso ~ 20M while the 
observer is set at r = 40M and we also assume that /3o ~ const. 




FIG. 10: The total energy emitted in gravitational waves {I = 2) from the homogeneous dust collapse with magnetic field for 
four values of initial radius Vso = 8M, 12M, 16M, and 20M as function of the magnetic field strength. 



C. Numerical results for the initial magnetic field profile (II) 

The second magnetic field profile considered in this work is the one that has its maximum close to the stellar surface 
and is described by equation (|116p . Here, again we assumed an initial profile for the fluid velocity that has Pq = const. 
and the initial radius of the dust spher is set to r^o = 8M . In Fig. 1111 we show the two components of the waveforms 



21 



^^^^ and for profiles (I) and (II). From this figure, we can see that the phase of l>(^) is almost independent of 
the initial profile of the magnetic field while the amplitude of the gravitational wave associated with the profile (II) 
is smaller than the one associated with the profile (I). 
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FIG. 11: Waveforms of quadrupolar gravitational radiation from the magnetized OS collapse associated with and "i>'^'', 
for the two initial profiles (I and II) for the magnetic field. 



In Fig. 1121 we compare the total energy emitted in gravitational waves as functions of the magnetic filed strength 
for the two initial magnetic field profiles. In this figure, we observe that the critical value etc, depends strongly on the 
initial profile of the magnetic field. Actually, for the profile (II) we get eI?^^^/{2M) = 3.267 x 10"^, Cgw.2/(2A/) = 
2.314 X 10~^, and etc = —7.08. These difference suggest that the form of the magnetic field affects in a critical way 
the amount of the emitted gravitational waves and worths more elaborate study. 
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FIG. 12: Total energy emitted in quadrupole gravitational waves from the homogeneous dust collapse with two different initial 
magnetic field profiles as functions of the magnetic field strength. The filled circles correspond to results for the profile (I) and 
the filled triangles for profile (II) 
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VII. CONCLUSION 

In this paper, we have studied the influence of magnetic fields on the eSiciency of collapse in emitting gravitational 
waves. Wc have considered an interior Oppcnhcimcr-Snydcr solution describing collapsing dust and wc studied how 
the amplitude and the waveform of the quadrupole axial perturbations were affected by the magnetic field, which 
actually enters as a second order perturbation term. These second order terms coming from the magnetic field are 
initially small but as the collapse proceeds they get amplified and become significant. For this study we have assumed 
that the magnetic field is axially symmetric and the Im = ^ magnetic field perturbations are the ones that couple 
to the I = 2 axial perturbations of the gravitational field. Additional assumptions have been made concerning the 
initial data and the influence of the magnetic fleld in the exterior. That is, we assumed momentarily static initial 
data independent of the magnetic field and we have not taken into account the influence of the exterior magnetic field 
in the propagating gravitational waves. 

The main result of this study is the proof of the strong influence of the magnetic field in the gravitational wave 
luminosity during the collapse. Depending on the initial profile of the magnetic field and its strength the energy 
outcome can be easily up to one order higher than what we get from the non-magnetized collapse. Additionally, we 
observed that the initial profile of the magnetic field perturbations can affect the energy output while it is possible to 
observe important phase shifts induced by the presence of the magnetic field. Since for a large initial radius the time 
needed for the black hole formation in longer, then the magnetic field acts for longer time on the collapsing fluid and 
its effect becomes more significant in the emitted gravitational wave signal. 

Concluding, we believe that although this study might be considered as a "toy problem" it has most of the ingredients 
needed in emphasizing the importance of the magnetic fields in the study of the gravitational wave output during the 
collapse. It is obvious that 3D numerical MHD codes will provide the final answer to the questions raised here, but 
this work provides hints and raises issues that need to be studied. 
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APPENDIX A: PERTURBED ENERGY-MOMENTUM TENSOR FOR THE ELECTROMAGNETIC 

FIELD INSIDE THE STAR 



The non-zero components of the perturbed energy-momentum tensor JT^i^*^'' for the electromagnetic field associated 



with dipole (Jm — 1) perturbations inside the star are given by the following relations 
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Following (29|, we can expand the perturbed energy- momentum tensor for polar parity perturbations in terms of 
tensor spherical harmonics as 
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where = Ptab + Hi + ^)Pnab/'2- [See, for the axial parity perturbations, equation (fT3|).] Then, it has been found 
that the non-zero tensor-harmonic expansion coefficients of ST^^^^^^ for Im — 1 are coupled with the I = and I — 2 
perturbations. 

The expansion coefficients for the I = perturbations are the following: 
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where (IA12[) belongs to axial parity perturbations, and (jA13|) and (|A14p belong to polar parity perturbations. 
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APPENDIX B: PERTURBED ENERGY-MOMENTUM TENSOR FOR THE ELECTROMAGNETIC 

FIELD OUTSIDE THE STAR 

The non-zero components of the perturbed energy-momentum tensor JT^i^*^' for the electromagnetic field associated 
with dipole {Im = 1) perturbations outside the star are given by the following relations 
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As for the interior, the non-zero expansion coefficients of STjxv'^^ for Zm = 1 electromagnetic field perturbations are 
coupled with the I = and 1 = 2 perturbations. 

The expansion coefficients for the I = perturbations are: 
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while the coefficients for the I = 2 polar parity perturbations have the form: 
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Finally, the coefficients for the I ~ 2 axial parity perturbations have the form: 
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